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Reading  on  this  talk 


•  www.stat.rice.edu/~Tiedi 

•  This  talk 

•  Intro  for  the  "untouched  mind" 

-  Explicit  computations  on  Binomial 

•  Monograph  on  "Multifractal  processes" 

-  Multifractal  formalism  (proofs,  references) 

-  Multifractal  subordination  (warping) 

•  Papers,  links 


Rudolf  Riedl  Pice  University 
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Why  Cascades 


Turbulence:  models  wanted 

•  Kolmogorov  1941  : 

<  [v(x+r)-v(x)p  >  ~  r^3 

•  Kolmogorov  1962  : 

<  [v(x+r)-v(x)p  >  ~  rH^ 

•  ...and  beyond 


Courtesy  Pr  Chainais 
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Measured  Data 


•  Networks 


•  Geophysics 


•  WWW 


•  Stock  Markets 


Multifractal  Analysis 


Toy  Example 
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The  Toy:  Binomial  Cascade 


•  Start  with  unit  mass 

•  Redistribute  uniformly 

portion  p  <Vi  to  the  left 
portion  1-p  to  the  right 


l 


0  ~2  1* 


Mass  re-distribution 


•  Iterate 

•  Converges  to  measure  jx 

oo 

t=  £fc/2A'  With  €fc  =  0, 1 

k=  1 

/(«i  ■••««)  :=  IE  E  a-/2fc  +  1/2")) 

k=  1  /e=l 

n 

ln(t)  :=  #{fc  <  n  :  ek  =  1}  =  JZ  ek 

k=  l 


p(/(«l  .  .  .  £n))  =  -  p)'”(t) 
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Multifractal  Spectrum 


Oscillate  ~  |t|a -»  local  strength  a  a=.l  a=.9  a=  g 


a(t)  lim  inf  an(t) 

n 


^  vJL  ^ 


Llliii  1.  llj.ii.iiL 


an(t )  = 


log  A /„(£) 


log  |/n(f)| 


I n ( 0 •  dyadic  interval  containing  t 


AIn(t ):  oscillation  indicator 
total  increment  over  /„, 
max  increment  in  In, 
wavelet  coefficients,... 


Collect  points  t  with  same  a  : 

r-.  r,  /,\  i  Dim(E,) 

Ea  —  {t  a(f)  =  a}  a 


Dim(Ea):  Spectrum 
-»prelevance  of  a 


.  _  I  £  r-s :  _  _l :  . 


_  i_  _ _ l.  .  / 


Binomial 


a(t)  :=  lim  inf  an(t) 


<*n(0  = 


log  A/ri(0 


We  take  dyadic  partition: 


log  |/n(OI 


ln(t)  —  #{fc  <  n  :  ek  =  1} 


«t)  =  iu l ...«»)  :=(E  E  «it/2*+l/2n)) 

fc=l  fc=l 


—  p^n(^)(l  _ p^n—ln(t) 


<*n(t)  =  -- — log2(p)  -  log2(l  -  p) 

n  n 

Range  of  exponents: 

t  0.  In  0,  ftn  *  Ooo  •—  IOg2(p)  ^  1 

t  =  1:  ln  =  n,  an  —>  a- oo  •=  -  log2(l  —  p)  >  1 
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"Typical"  exponents 


Recall 


t=0,  t=l  seem  "atypical". 
Intuition:  for  a  "typical"  t: 


a(t)  :=  lim  inf  an(t) 

n 


an(t)  = 


log  A/„(Q 

log|/n(OI 


ln(t)  —  #{fc  <  n  :  €fc  =  1} 


ln(t)  ~  n/2 

Rigorously:  Law  of  Large  Numbers 

•  Binary  digits  ek  are  independent,  P[ek=0]=  P[ek=l]  =  Vi\ 

•  t  is  uniformly  distributed  (i.e.,  with  Lebesgue  measure  Z) 


ln(t) _ ^  ip  r  i _ i/o 

n  nk=  1 


"Typical"  exponent 


<On(0  — 


- — log2(p)  -  log2(i  -  p) 


n  n 

a0  := --|og2(p)  - -!°g2(l -p)  >1 
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A  first  point  on  the  Spectrum 


Conclusion: 

*  ^  ^ 

•  Mass  Distribution  Principle 

(Lebesgue  measured  is  1-dim  Hausdorff  measure) 


dim£'an  =  1 


"Where"  and  "how  many"  are  the  other  exponents? 
•  Choose  digits  "unfairly",  e.g.,  prefer  1  over  0. 
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Other  exponents 

The  measure  jll  prefers  1  over  0  (ratio  1-p  to  p). 
Intuitive: 

ln(t)  ~  n(  1  -  p) 

Rigorously:  Law  of  Large  Numbers  using  p 

•  Binary  digits  e  are  independent,  P[ek=0]=p,  P[ek=l]  =  1-p: 

•  t  is  distributed  according  to  p 


Recall 


etn(0  = 


log  A/n(0 


log  \In(t)\ 

ln(t)  =  #{fc  <  n  :  ek  =  1} 


n 


1  n 

=  -  2J  ek  -*■  E4el  =  1  -  P 
71  k= 1 


•  p-typical  exponent _ 

,  x  n  ~  Tn^T)  ~  "  /n(0  .  T  T 

an(t)  = - iog2(p) - log2(l  -  p) 

n  n 

_ ->  ai  := -plQ92(p)  -  (1  ~P)  <092(1  ~P)  <1 
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A  second  point  on  the  Spectrum 


Conclusion: 

•  >  0 


(Hausdorff  dimension  of  jli?  It  is  a1<l!) 


log //(/„(/)) 

M°  =  iog]7^(0]  * 01 


Large  Deviations 

and  the 

Multifractal  Formalism 
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Box  Spectrum 


Recall 


Notation: 


/ri(0  —  I  (^1  •  •  •  ^n) 

log  A In(t) 


log  |/n(Ol 

Eg  :=  {t  :  a(t)  =  a} 


Nn(a ,  (5)  :  =  #{(ei . . .  en)  :  a— 8  <  an(e \  . . . en)  <  a+5}. 


/(a)  :  =  lim  lim  sup  -  log2  iVn(a,  8) 
<5|0  n—y  oo  n 


Proof 


Thm:  we  always  have 


Fix  n  To  Drove  the  first  Dart  of  the  lemma 
consider  an  arbitrary  t  >  /(a),  and  choose  n  > 
0  such  that  i  >  /(«)  +  2t),  Then,  there  is  e  >  0 
and  integer  m0  such  that 

.V„(u,r)  < 

for  all  >■  >  triQ.  Let  us  define  J(m)  :=  u(L„  : 
n  >  ni  and  a  -  c  <  n*  <  o  +  «}.  Then,  for  any 
m  the  intervals  /J  with  kn  €  J(rn)  form  a  cover 
of  £>.  These  Intervals  are  of  length  less  than 
2  Moreover,  for  any  m  >  mo  we  have 


dim Ea  <  /(a ) 


E  IGP  =  E  *»(«.*>  2-« 

^  <  £  2-«V 

n>f*»  H>lll 

tends  to  zero  with  m  —  oo.  We  conclude  that 
the  -.-dimensional  Mausdorff  measure  of  Ea  is 
zero,  hence.  dlmJL'„  <  Letting  7  — .  /(«) 
completes  the  proof 


[www.stat.rice.edu/~riedi] 


Beware  the  folklore:  f(a)  is  NOT  the  box-dim  of  E 
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Recall 


Legendre  spectrum 


Notation:  partition  sum  and  function 


Ai(0  —  I  (^1  •  •  •  ^n) 
log  A In(t) 


C*n(  0  — 


log  |/n(0l 

Eg  :=  {t  :  a(t)  =  a} 


Sn(q):=  Y.  |A/n(ei...en)|9=  £  |2n|9an(ei-en) 

f  1  ^1 


r(g)  :=  lim  inf  —  log2  Sn(q) 

TL - >00  rjr^ 


•  Thm:  we  always  have 

/(a)  <  r*(a)  :=  inf(ga  -  t(<?)) 

Proof 

£  |A/„(£l...€„)|«  >  £  |A/„(n...e„)|«  >  iV„(a15)2-n(<?a+W<5) 

Cl"*Cn  an(ci*..Cn)€[a-5,a+<5] 

>  2~n(9a-/(a)+<5'+M<5) 


Rudolf  Riedi  Rice  University 


stat.rice.edu/~riedi 


Legendre  spectrum 

•  Thm:  provided  an(t)  are  bounded  we  have 

/(a)  =  r*(a)  for  a  = 


•  Proof  idea:  steepest  ascent  (large  deviations) 

m 

'V*.  |A/n(ei . . .  en)|f/  <  y:  ^  ,  |A/n(€i...en)|9 

Cl—Cn)  *=1  a„(e1...c„)G[/<5-<5,Z<5+<5] 

m 

<  Y.  Nn(lS,S)2~n^16-^ 

1=1 

<  jr  2-n(^'-/W-<5  ~M<5)  <  m2-n(info(ga-/(a)-<yi  -|q|<5) 

/=1 

-Thus:  t(q)  = /*(q)  =  inf(ga  - /(a))  for  all  q. 

-  r  is  concave,  non-decreasing,  differentiable  with  exceptions 

-  Recover  f=f**  at  a=r'(q)  using  lower  semi-continuity 
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T(q)-> 


Legendre  transform  101 

•  Elementary  calculus: 

r*(a)  :=  inf (qa  -  r(q))  =  qa  -  r(q) 
where  q  is  defined  by  a  =  r'(q ) 

-  Tangent  of  slope  a  to  r(q) 

-  Intersection  with  ordinate  yields  -7*(a) 


Binomial  Spectrum 


continued 
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Partition  function  of  the  Binomial 


Sn(q)  —  | ^ -^n(ci  •  •  •  €n) 

=  S  lPn~ln  (ei...en)  (1  —  pY 

=  b«  +  (i-P)«]". 

•  (Upper)  envelop  of  dim(Ea): 

i 

r(q)  =  liminf  —  log2S,n(g) 

71 — ►  OO 

=  -  log2 [p9  +  (l-p)f/] 


n(ci...€n) 
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Insight  from  Large  Deviations 


•  From  steepest  ascent: 


-n(\r\fu(qa-f(a)) 


Sn(q)  =  £  |A/n(ei...en)|9^2 

C 1  •• 

=  2-n^a-/(«))~  £  lA/nCei-.-en)!9 

ttn  Rl  ® 

Dominant  terms  in  Sn(q),  for  fixed  q,  are  the  ones  with 

log  A/n  _  , 

an(ei  . . .  en)  =  - - —  ~  a  =  r  (q) 

log  \In\ 

•  ...and  vice  versa:  these  terms  contribute  such  that 


Sn(q)~  2~nr^  =  (p«+(l-p)9) 


n 


For  the  Binomial  these  correspond 
to  mass  re-distribution  in  ratio  pq  to  (l-p)q 
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Locating  the  exponents 


Fix  q. 

Consider  the  measure  j u,q  defined  as  p.  but  with  mass  ratio  to 
(l-p)q.  Intuitively,  we  have  then: 

ln(t)  =  #{fc  <  n  :  ek  =  1}  ~  n  rz  =  n(l-p)<?2r(<?) 

pq  +  (1  -  V)q 

Rigorously:  Law  of  Large  Numbers  using  jLiq 

•  Binary  digits  e:  indep,  P[ek=0]=pq2x(q),  P[ek=l]  =  (l-p)q2x(q) 

•  t  is  distributed  according  to  juq 

n  n ,  , 


•  jLt-typical  exponent 

acn(t)  -  -n-lnty  iog2(p)  -  ^^log2(l  -p) 

n  n 

_ — »  aq  :=  -pq2T^  \OQ2 (p)  —  (1  -  p)q2T^  log2(l  —  p) 
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Completing  the  Spectrum 


Conclusion: 

•  f-LqiEaq)  >  o 

aq  =  r\q ) 


Recall 

r(q)  =  ~  \og2[p(l  +  (1  -  p)f/] 

«„  =  -p,'2r(,»>log2(p)-(l-p),'2T(")log2(l-p] 


•  Hausdorff  dimension  of  |j,q : 


lOQ 

log  |/nCOI 


n  1 \og2[pq2T(q)] 

n 


=  -r(q)+qan(t) 

->  qaq  -  r(q)  =  T*(aq) 


—  log2[(l-p)92T^] 

n 


•  Mass  Distribution  Principle 


Rudolf  Riedi  Rice  University 


dimgga  >  r*(qg) 
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T(q)-» 


Lessons 


Binomial  cascade:  6\rr\Ea  =  /(a)  =  r*(a) 


•  Points  with  exponent  logju(I(t))/log|I(t)|  ~  a=r'(q) 

-  Are  concentrated  on  the  support  of  j uq 

-  Dominate  the  partition  sum  Sn(q) 

•  Partition  function  allows  to  bound/estimate  dimfE.,) 
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Random  Cascades 


A  further  multifractal  envelop 
Convergence  and  Degeneracy 
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Multifractal  Spectra  and  Randomness 


•  Aln(t):  oscillation  indicator  for  process  or  measure 

•  Pathwise  Keca" 


dimBa  <  /(a)  <  r*(a) 


Ea  —  {t  :  a(f)  =  a} 
Nn(a,e)  ~  2n^a'> 


•  Sn(q)  is  q-th  moment  estimator. 

•  Replace  by  true  moment: 

T(q )  :=  lim  inf-- log2E5n(g) 

71 — ►  OO 


SnU l)  ^  2-"T^> 

Recall 

Sn{(j)  =  ^  l^4i(€l  •  •  •  en)\q 

_ gj...  fn _ 


•  ...analytically  easier  to  handle  and  often  sufficient 

•  T(q)  is  concave  like  r(q),  but  NOT  always  increasing 
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Pathwise  and  deterministic  envelop 

•  Lemma:  With  probability  one  for  all  q  with  T(q)  <  oo. 

r(q,uj)  >  T(q ) 

[Proof:  www.stat.rice.edu/~riedi] 

•  Cor:  r*(a, cu)  <  T*(a ) 

E|log(X)|  <  log 

•  Weaker  result  from  Chebichev  inequality: 

E[r(q,u;)]  >  T(q ) 

Quenched  Average  Annealed  Average 

•  Material  science:  free  energy  is  "self-averaging"  iff 
quenched  and  annealed  averages  are  equal. 
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Multifractal  Envelops 


•  Almost  surely,  for  all  a: 


Recall  at  a  =  r'(q ) 

/(a)  =  t*(q) 


dim£a  <  /(a)  <  r*(a)  <  T*(a) 


•  Holds  always  provided  use  same  AIn  in  all  spectra 

•  Choice  of  scales  In 

-  In  is  here  dyadic,  could  be  any  sub-exponential 

-  This  could  affect/change  f,  r  and/or  T  due  to  boundary  effects 

-  Robust:  Aln  =  oscillation  in  In  and  its  neighbor  intervals 

•  Choice  of  oscillation  indicator  Aln 

-  For  true  Hoelder  regularity  Aln  =  max  increment  "around"  In 

-  Aln  =  Wavelet  coefficient:  only  a  proxy  to  Hoelder  regularity! 

-  For  measures  supported  on  [0,1]:  A ln  —  n(ln)  gives  Hoelder! 
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Multifractal  Envelops 


•  Almost  surely,  for  all  a: 


Recall  at  a  =  r'(q) 
/(a )  =  T*(a) 


dim£a  <  /(a)  <  r*(a)  <  T*(a) 


•  Special  feature: 

-  If  a  property  of 
"bounded  total  variation" 
holds  then  the  spectrum  f 
touches  the  bi-sector: 


It  Eei...e„  AJn(ei . . .  e„)  <  C  for  all  n 
then  t(1)  =  0. 
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Multifractal  Envelops 


•  Almost  surely,  for  all  a: 


Recall  at  a  =  r'(q ) 

/(a)  =  t*(q) 


dim£a  <  /(a)  <  r*(a)  <  T*(a) 

•  Terminology: 

-  Multifractal  formalism  "holds"  if 

•  dim(Ea)=f(a)=T*(a) 

with  your  preferred  oscillation  indicators  Aln, 
e.g.,  Holder  exponent  in  Ea,  wavelet  decay  in  f(a). 

[First  step:  show  T  is  same  for  Holder  and  wavelets.] 

-  Falconcer:  "A  concise  definition  of  a  multifractal  tends  to  be 

avoided." 

-  Others:  "An  object  is  multifractal  if  the  formalism  holds  for  it." 

-  Others:  "An  object  is  multifractal  if  it  has  more  than  one 
singularity  exponent",  (not  mono-fractal) 
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Multifractals 
and  classical  regularity 
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Besov  spaces 


•  For  oscillation  indicator  from  wavelets: 

sup {s  :  Y  e  BSV{LU )}  =  1 

u 

•  Proof:  use  wavelet  coefficients  Cj  k=  Alj(k2J)  and 
equivalent  Besov  norm 

/  \  l/v  /  \  v/u\  1/*’ 

El  A), or  +  E  ( E  2jsu2-i  \2^2Cjtk\v  ) 

V  k  )  \j>J o  V  &  J  } 
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T(q)-> 


Kolmogorov 

•  Thm  [Kolmogorov]: 

—  If  E[  |  A(s)-A(t)  |b  ]  <  C  |  s-t  |1+d  then  almost  all  paths  of  A 
are  of  (global)  Holder-continuity  for  all  h  <  d/b, 

-  i.e.,  for  all  h  <  T(q)/q  . 

•  The  best  such  h  is  min(a  :  T*(a)>0). 

-  T(q)/q  =  slope  of  tangent  through  the  origin. 


F  q  -> 


Binomial  Spectrum 


continued 
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Binomial  with  Random  Multipliers 


•  Random  re-distribution  i - 

•  Multipliers  Independent  | _ * 

between  scales  0  V2 


/J-rz ( -^ (€  1  •  •  •  en))  —  . 


•  Conservative: 

Me\...enO  “1“  ^ei...enl 


—  1  Conservative  re-distribution 


•  Conservation  is  too  restrictive 
for  stationarity! 

•  "Martingale  de  Mandelbrot": 

^[■^ei...enO  ”1“  ■^ei...€nl]  —  ^ 


1 


0  ~2  1* 


Independent  re-distribution 

stat.rice.edu/~riedi 
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Convergence  of  Random  Binomial 


•  Conservative: 

M ei-..enO  d-  ^e\...€nX 

-  For  all  m>n 

1  •  •  •  en))  —  -Wei  •  •  • 

-  Thus  converges  to 

•  •  •  en))  —  A/gj  .  .  .  Me\...en 
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Convergence  of  Random  Binomial 


•  "Martingale  de  Mandelbrot": 


-  A  price  to  pay  towards  stationarity  « 


-  Martingale:  For  all  m>n 

(ei  . .  .  €n)) l-T'n]  —  =  (el  •  •  •  en)) 

-  Thus  converges  almost  surely  (but  may  degenerate) 

-  We  have 


[yrt( V (ei  . . .  £n))|-?"" rz]  —  •  ■  •  -Wei 
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Envelope  for  Random  Binomial 

•  By  independence  of  multipliers 

-  Martingale  of  Mandelbrot: 

E[Sn(<7)]  =  Z  E|A/n(ei...e»)|9=  Z  E|Mei . . .  M£l...en|9  =  2nE[M9]n. 


T(q)  =  -1  -  log2E[M9] 

-  Conservative:  similar 

E[Sn(<?)]  =  Z  H[M9]n_Wei-€n)H[(l  -  M)g]<n(€l"  €n) 

=  (TE[Mq]  +  E[(l  —  M)9])n 
-  (2E[M9])n. 
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Kahane-Peyriere  theory  for  the 
Martingale  of  Mandelbrot 


Martingale  "degenerates" 

-  iff  |j,([0,l])=0  almost  surely  zero 

-iff  E  |ii([0,l])=0 

-  iff  T'(1)<=0 

Intuition: 


-0.5  0  0.5  1  1.5 


-  T'(l)  =  =  dimension  of  the  carrier  of  jli. 

-  If  T'(1)>0  then 


•  3q>l  with  T(q)>0 

•  ju.  converges  in  Lq 

•  E[/x([0,l])]  =  lim E[/xn([0, 1])]  =  1 
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Multifractal  formalism  holds 


•  Thm  for  random  binomial  [Barral,  Arbeiter-Patschke,  Falconer]! 

-  Set  AIn  =  ju(In). 

-  Assume  M  has  a  finite  moment  of  some  negative  order 

-  Then,  with  probability  1:  for  all  a  such  that  T*(a)>0 


dim£a  =  /(a)  =  r*(a)  =  T*(a) 


•  Note: 

-  T*(a)>0  means  a=T'(q)  with  q  limited  by  tangents  through 
the  origin:  T'(q)=T(q)/q. 

-  Little  known  in  general  for  other  a  ...or  q!  Possible:  r(q)>T(q) 

-  Proofs:  Use  Mass  distortion  Principle  with  factors  Mq 
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Wavelets  for  the  Binomial 


•  Compactly  supported  wavelet 

-  AIn  =wavelet  coefficient  corresponding  to  In 

-  AIn  same  rescaling  property  as  measure  itself 

-  Same  T(q) 

-  Multifractal  formalism  holds 
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Toy  examples 
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Log-Normal  Binomial 


Deterministic  envelope  is  a  parabola:  [Mandelbrot] 


T(g)  =  (g- 1)  1  - 


0" 


2  ln(2) 


g)  for  q  <  gcrit  :=  2  ln(2)/cr2 


•  Zeros:  q  =  l,  q=qcrit 

•  Non-Degeneracy:  .  0 

T  (1)  >  0  ^  gcrjt  >1^2  I  n  (2)  >  cr2 


•  Spectrum  is  parabola  as  well 


T*(a) 


i  _  lnC2)  /  <t2  \2 

2cr2  v'  2  ln(2) J 


•  Partition  function  x(q)  is  non-decreasing, 

•  thus  i(q)  >  T(q)  (at  least)  for  q>(l+qcrit;)/2 
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Multifractal  Product  of  Pulses 

together  with 

I.  Norros  and  P.  Mannersalo 
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Network  Traffic  is  Multifractal 


•  Visually  striking 

•  Scaling  of  impressive  quality 

(Levy  Vehel  &  RR  '96, 

Norros  &  Mannersalo  '97, 

Willinger  etal  '98) 

•  Statistical  models: 

-  Binomial  cascades  with 
scale  dependent  multipliers 

(Crouse  &  RR  '98,  Willinger  et  al  '98) 

•  Not  stationary! 

-  Cumbersome  for  statistics 

-  and  probability  (Queueing) 


tsooo 
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Multifractal  paradigm 

Multiplicative  Processes: 

-  From  redistributing  mass  to  multiplying  pulses 

A(t)=  lim  I  An(s) ...  An(s)ds 

n— »oo  7q 


Binomial  Cascade 

-  An(s)  is  constant  on  dyadic 
intervals 

-  Conservative: 

An(2k/2")  +  An((2k+l)/2")  =2 

-  Martingale  de  Mandelbrot: 

E  An(s)  =  1 

-  Not  stationary 


Ai(s) 


i  2M  - 

1 . . ► 


_  t 

. . ► 

0  1/4  1/2 
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Multifractal  paradigm 

•  Multiplicative  Processes: 

A(t)=  lim  I  An (s)  ...  An(s)ds 

71— >  OO  J  Q 

4 

it - 


•  Stationary  Cascade 

■ 

-  An(s)  is  stationary  If 

-  Conservation: 

EAn(t)  =1  f 

-  "self-similarity":  1t 

An(s)  =d  A^sb11)  .1. 
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Parameters  and  Scaling 

•  Parameter  estimation 

-  Aj(s):  i.i.d.  values  with  Poisson  arrivals  (Aj): 

-  Z(s)  =  log  [  Ai(s)  A2(s)...  An(s)  ] 

-  Cov(Z(t)Z(t+s))=  Z)j=1  n  exp(-A|S)Var  A|(s) 

•  Performance  of  predictors  /  simulations 


0.005  . 


A-l  A2 


0.00<( 


Multifractal  Envelope 

(with  Norros  and  Mannersalo) 

T(q)=q-l-log2E[Atl] 
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Interlude 


Self-similar  processes 
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Statistical  Self-similarity 

•  Self-similarity:  canonical  form 

-  B(at)  =fdd  C(a)  B(t)  B:  process,  C:  scale  function 

-  Iterate:  B(abt)  =fdd  C(a)C(b)  B(t) 

-  C(a)C(b)=C(ab) 

->  C(a)  =  aH  :  Powerlaw  is  default 

•  H-self-similar: 

B(at)  =fdd  aH  B(t) 


•  Examples 

-  Gaussian:  fractional  Brownian  motion  BH(t)  is  unique  H-self- 
similar  Gaussian  process  with  stationary  increments. 

-  Stable:  not  unique  in  general,  a=l/H:  Levy  motion 
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Statistical  Self-similarity 

•  How  do  self-similar  processes  occur? 

-  Xk:  stationary  time  series 

-  U(t):=X1  +  ...  +  X[t] 

-  If  U(nt)/f(n)  f  d  d  Z(t) 

-  then  necessarily  H=limn^1D  log  f(n)  /  log(n)  exists 
and  Z(t)  is  H-self-similar. 

-  If  Xk  are  iid  with  finite  variance, 
then  H  =  l/2  and  Z  is  Brownian  motion 

-  If  Xk  are  LRD,  then  H>l/2  and  Z  is  fractional  Brownian  motion 

•  Prediction  and  estimation  windows 
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Self-similar  Processes 


•  What  do  they  model? 

-  Sustained  excursions  above/below  the  mean 

•  Different  from  (finite  order)  linear  models 

-  Auto- Regressive 

-  ARMA 

-  (G)ARCH 

-  Exponential  decay  of  correlations 


Corresponds  to  infinite  order  AR  models 


-  FARIMA 

-  FIGARCH 
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fBm(t)  =  _J  K(t,s)  dW(s) 
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Multifractal  Subordination 

Processes  with 
multifractal  oscillations 
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Multifractal  time  warp 

BH(M(t)):  Bh  fBm,  dM  independent  measure 


A  versatile  model 

-  M(t):  Multifractal 
Time  change 
Trading  time 

-  B:  Brownian  motion 

Gaussian  fluctuations 
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Holder  regularity 


•  Levy  modulus  of  continuity: 

-  With  probability  one  for  all  t 

\BH(t  +  S)-BH(t)\~\S\H 


-  Thus,  exponent  gets  stretched: 

\BH(M(t+6))-BH(M(t))\  ~  \M(t+5)-M(t)\H  ~  \S\Ha^ 


Multifractal  formalism  for  BH(M(t)) 


•  Conditioning  on  M  one  finds: 

~  |<j| 

'thUS  TB(M)(q)  =  TM(qH) 

-  which  confirms  the  stretched  exponent: 

TB(A ofo)  =  ht'm^H) 

-  and  matches  with  warp  formula  before: 

tb(m)^  =  t Ua/H) 

-  If  the  formalism  holds  for  M,  then  also  for  BH(M(t)) 
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Auto-Correlation 

•  Conditioned  on  knowing  M(t): 

-  E[  B(M(t))  B(M(s))  I  M  ] 

=  (a2/2)  [  M2H(t)  +  M2H(s)  -  M2H(t-s)  ] 

-  Non  stationary  Gaussian  Process 

-  Increments:  X(t)  =  B(M(t+l))  -  B(M(t)) 

-  E[  X(t)  X(s)  |  M  ]  =  (o2  /2)  x 

(  [M(t+1)  -  M(s)  ]2H  -  [  M(t)  -  M(s)  ]2H  -  [M(t+1)  -  M(s+1)  ]2H  +  [  M(t)  -  M(s+1)  ]2H ) 

•  Unconditioned:  For  H  =  l/2  and  E[M(t)]=t 

-  E[  B(M(t))  B(M(s))  ]  =  a2  min(s,t) 

-  E[  X(t+k)  X(t)  ]  =  E[  M(k+1)  -  2M(k)  +  M(k-l)  ]  =  0 

-  Uncorrelated  increments,  stationary,  2nd  order 

-  But  not  Gaussian 

-  Dependence  of  higher  order  through  M(t) 
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Estimation:  Wavelets  decorrelate 


(with  P.  Goncalves) 

•  wjk  =  J  v|/jk(t)  B(M(t))  dt 

N:  number  of  vanishing  moments 

*  E[WJk  WJm  ] 

=  J  J  'PjKtt)  f)m(s)  E[B(M(t))  B(M(s))]  dt  ds 

=  J  J  V‘)  'J'jm(s)  E[|M(t)  -  M(s)|2H]  dt  ds 

~  0(  |k-m|T<2H>  +1  -  2N  )  (  |k-m|->  °o) 
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Multifractal  Estimation  for  B(M(t)) 

•  Weak  Correlations  of  Wavelet-Coefficients: 

(with  P.  Goncalves) 


08 
0  6 
04 
02 
0 

-02 
-0 

•  Improved  estimator  due  to  weak  correlations 

•  Multifractal  Spectrum 

M(t+s)  -  M(t)  ~  sa^) - 

B(t+u)  -  B(t)  ~  uH  (v  t) 

B(M(t+s))  -  B(M(t))  ~  sH*aW 


Daubechies2 
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Estimation 
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From 

Multiplicative  Cascades 

to 

Infinitely  Divisible  Cascades 


with 

P.  Chainais  and  P.  Abry 
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Independent  work: 
Castaing,  Schmidt, 
Barral-Mandelbrot,  Bacry-Muzy 
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tog2(Variance) 


Adapting  to  the  real  world 

Real  world  data 

•  can  deviate  from  powerlaws:  traffic 

•  has  no  preference  for  dyadic  scales 

Lukacs:  if  the  data  does  not  fit  to  the 
model  then  too  bad  for  the  data. 
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Experimental  results 

Courtesy  P.  Chainais 


Irt  Conf  96% 


H(q)  n(a):  non-powerlaw 
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Beyond  Self-similarity 

•  Self-similarity  revisited: 

-  B(at)  =d  C(a)  B(t)  B:  process,  C:  scale  function 

-  B(abt)  =d  C(a)C(b)  B(t) 

-  C(a)C(b)=C(ab)  ^  C(a)  =  aH 

-  E [  |  B(an)  1^]  =  c(q) 

-  linear  in  q  (mono-fractal) 
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Beyond  Self-similarity 

•  Self-similarity  revisited: 

-  B(at)  =d  C(a)  B(t)  B:  process,  C:  scale  function 

-  B(abt)  =d  C(a)C(b)  B(t) 

-  C(a)C(b)=C(ab)  ^  C(a)  =  aH 

-  E [ |  B(an)  I0)]  =  c(q)  (aciH)n 

-  linear  in  q  (mono-fractal) 

•  More  flexible  rescaling  "Ansatz": 

-  C=C(a,t)  ?  :  non-stationary  increments 

-  C=independent  r.v.  for  every  re-scaling  : 

-  X(a...at)=  X(ant)  =  C1(a)...Cn(a)  X(t):  multiplicative 
-E[|  X(an)  |q]  =  c(q)  E[|  C(a) 

-  non-linear  in  q;  powerlaw 
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Infinitely  divisible  scaling 


Self-similarity:  IE[|.B(t  +  5)  —  B(t)\q]  ~  5qlf 
Multifractal  scaling:  TE[\M  (t  +  5)  —  M  (t)\q]  ~  S1+T^ 

IDC  scaling:  E[|X(t  +  S)  -  X(t)\q]  cs  exp[n(<5)C(g)] 


•  Multifractal  scaling  reduces  to  self-similarity  if  T  is 
linear  in  q.  (sometimes  called  mono-fractal) 

•  IDC  reduces  to  multifractal  scaling  if  n(S)=-log(6) 

•  In  general  n(S)  gives  the  speed  of  the  cascade 
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Geometry  of  Binomial  Pulses 


Time-Scale  plane:  codes 
shape  of  pulses 

-  Position  (T=center) 

-  Size  (R=length) 
Pulses: 

Pj(t)  =  Wj  if  |t-t,|<r,/2 

1  else 
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For  Binomial: 
Strict  dyadic 
geometry 
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Stationary  geometry 


Large  Scale  Pulses 


Medium  Scale  Pulses 
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Compound  Poisson  Cascade 


Poisson  points  (tj  ,  r;  )  in  time-scale  plane  with  marks  Ws 
Cone  of  influence  at  t 

C(r,t)  =  {(ti,rj)  \t-ri<ti<t,ri>r} 

Cascade  Process: 

Qr(t)  =  I]  Wi 

( ti,n)eC(r,t ) 

r 

•  Poisson  Cascades  exhibit  scaling  properties  akin  to  IDC  scaling 

m(C(r,t ))  =  m(C(r,  0))  =  E[#{Ri,rj)  €  C(r,t)}] 

IE Qr(t)q  =  exp  [-(p(q)m(C(r,  *))] 
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Cascade  and  AR  processes 


•  Continuous  version  (IDC): 

Q(t )  =  exp 

=  exp  j  k(j(t,  s)dM(s) 

-  M  is  an  infinitely  divisible  measure 


•  Classic  theory  to  be  exploited: 

-  AR-type  processes 

B„(t)  =  jkH{s,t)dW(s) 


-  kernel  estimate  of  the  random  measure  dM 
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Cascades:  Invariance  and  scaling 


Infinitely  divisible  nature 
and  scaling  of  the  cascade: 

Qr(t)  =  n  wi  =  Hwixl[wi 


C(r,t ) 


Cb 


cl 


Cb(t) 


=  QbWxHWi 

Cr 

Rescaled  version  of  Qr/b 
in  the  scale-invariant  case  only! 
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Poisson  Cascade  has  re-scaling  properties; 
in  scale  invariant  case:  akin  to  Product  of  Processes 
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Multifractal  scaling 


Multifractal  formalism  holds  in  self¬ 
similar  case  [Barral-Mandelbrot] 

Infinitely  Divisible  Scaling 


Recall 


IQrjty1  =  exp  [-tp(q)m(C(r,  *))] 


IEy4(f)9  ~  exp  [— *)) 

-  powerlaw  only  if  m(C(t,*))  =  -log(t) 

-  for  IDC  in  Self-Similar  Case  [Bacry-Muzy,Barral] 

-  for  CPC  and  log-normal  IDC  in  certain  non- 
powerlaw  cases  [Chainais-R-Abry] 
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Simulations 


•  Stationary  Cascade: 


•  Non-powerlaw  scaling 
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"Never  happy":  More  flexibility 


•  Better  control  of  scaling 

•  Wider  range  of  known  non-powerlaw  scaling 


O  © 

0  0,  .c 

•  Higher  dimensions:  anisotropy  V  V 

-  "As  expected"  in  generic  cases  [Falconer,  Olsen] 

-  Formalism  may  break  if  directional  preferences 


[McMullen,  Bedford,  Kingman,  R] 
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Overall  Lessons 

•  Multifractal  spectrum  <->  regularity 

-  Besov  spaces 

-  Global  Hoelder  regularity 

•  Powerful  modeling  via  multiplication  through  scales 

-  Poisson  product  of  Pulses 

-  Multifractal  warping 

-  Degeneracy:  price  to  pay  for  stationarity 

•  Estimation  via  wavelets 

-  Multifractal  envelopes 

•  numerical  x(q), 

•  Analytical  T(q) 

-  Choice  of  wavelet,  of  order  q 

-  Interpretation:  what  kind  of  spectrum  did  you  estimate 

•  Hoelder  exponent 

•  Wavelet  decay 
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To  take  away 


•  Cascades  matured  to 
versatile  multifractal  models 


•  There  remains  much  to  do. 
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Reading  on  this  talk 


•  www.stat.rice.edu/~riedi 

•  This  talk 

•  Intro  for  the  "untouched  mind" 

-  Explicit  computations  on  Binomial 

•  Monograph  on  "Multifractal  processes" 

-  Multifractal  formalism  (proofs) 

-  Multifractal  subordination  (warping) 

•  Papers,  links 
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